Center Vortices and Chiral Symmetry Breaking in SU(2) Lattice Gauge Theory 



o 

U 

Oh 
< 






> 

(N 

o 



X 

5-1 



Roman Hollwiesei^ Thomas Schweigler, and Manfried Fabeijf] 



Institute of Atomic and Subatomic Physics, Vienna University of Technology, Operngasse 9, 104-0 Vienna, Austria 

Urs M. Heller 

American Physical Society, One Research Road, Ridge, NY 11961, USA 

(Dated: April 5, 2013) 

We investigate the chiral properties of near-zero modes for thick classical center vortices. In par- 
ticular we analyze the creation of near-zero modes from would-be zero modes of various topological 
charge contributions from center vortices. We show that colorful spherical vortex and instanton 
ensembles have very similar Dirac eigenmodes and also vortex intersections are able to give rise 
to a finite density of near-zero modes, leading to chiral symmetry breaking via the Banks-Casher 
formula. We discuss the influence of the magnetic vortex fluxes on quarks and how center vortices 
may break chiral symmetry. 
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I. INTRODUCTION 

Quantum chromodynamics (QCD) at low energies is dominated by the nonperturbative phenomena of quark con- 
finement and spontaneous chiral symmetry breaking (%SB). Presently, a rigorous treatment of them is only possible 
in the lattice regularization. Many of the important features of non-Abelian gauge theories are already present in 
SU(2), which simplifies theoretical and numerical calculations. 

The nonperturbative vacuum can be characterized by various kinds of topological gauge field excitations. A well 
established theory of xSB relies on instantons [IH1] , which are localized in space-time and carry a topological charge 
of modulus 1. According to the Atiyah-Singer index theorem [SHE], a zero mode of the Dirac operator arises, which is 
concentrated at the instanton core. In the instanton liquid model [§1-111| overlapping would-be zero modes split into 
low-lying nonzero modes which create the chiral condensate. 

Center vortices [12H17J . on the other hand, are promising candidates for explaining confinement. They form closed 
magnetic flux tubes, whose flux is quantized, taking only values in the center of the gauge group. These properties 
are the key ingredients in the vortex model of confinement, which is theoretically appealing and was also confirmed by 
a multitude of numerical calculations, e.g. [TH1HS]- Lattice simulations indicate that vortices may be responsible for 
topological charge and xSB as well 20 25J , and thus unify all nonperturbative phenomena in a common framework. 
A similar picture to the instanton liquid model exists insofar as lumps of topological charge arise at the intersection 
and writhing points of vortices. The colorful, spherical SU(2) vortex as introduced in previous articles of our group 
[2l?l - f2T)] . may act as a prototype for this picture, as it contributes to the topological charge by its color structure, 
attracting a zero mode like an instanton. 

In this article we want to show how the interplay of various topological structures from center vortices (and 
instantons) leads to near-zero modes. Using the overlap and asqtad staggered Dirac operator, we compute a varying 
number of the lowest-lying Dirac eigenf unctions, including the zero modes. By visualizing the probability density, we 
compare the distribution of the eigenmode density with the position of the vortices and the topological charge created 
by intersection points and color structures. These results manifest the importance of center vortices also for chiral 
symmetry breaking. 

II. TOPOLOGICAL CHARGE FROM CENTER VORTICES 

From the definition 

Q OC tuvapFpvFap , (1) 

it is clear that topological charge contributions arise where two perpendicular nontrivial plaquettes meet. On a center 
vortex, the vortex sheet is locally orthogonal to the nontrivial plaquettes. Hence, topological charge emerges at so- 
called singular points of center vortices, where the set of tangent vectors to the vortex surface spans all four space-time 
directions |30j . There are two possibilities where this can occur: 

• intersection points of two different surface patches: Q = ±1/2 , 

• writhing points of a single surface patch: \Q\ < 1/2 . 

In lattice language, intersection points have four plaquettes attached to them extending in two particular space time 
directions (e.g. xy), and another four plaquettes extending in the two orthogonal directions (e.g. zt). By contrast, 
at a writhing point one can start at any attached plaquette and pass around over all others along a continuous path 
since they are connected by common links. 

In SU(2), a Wilson loop cannot distinguish between center fluxes of opposite sign, since exp(i-7r) = exp(— i7r) = — 1. 
However, the sign of the topological charge is sensitive to the direction of the flux. To preserve this information for 
thin vortices, one can assign an orientation to the vortex patches. Note that this quantity is not directly linked to the 
geometry of the vortex. A nonorientable vortex need not be nonorientable as a surface [31] . 

Vortices always have closed surfaces. Furthermore, lattice studies show that the major fraction of vortex patches 
belongs to a single large vortex winding through the whole lattice. Intersections between different vortices are therefore 
only of minor importance. But for a closed oriented surface the total charge is zero. Charged vortices must therefore 
be globally nonoriented and consist of differently oriented patches, which are separated by monopolc worldlines. The 
colorful spherical vortex is a classical representation of such a vortex configuration, contributing to topological charge 
only through its color structure. It will be discussed shortly in Sec. |IIB| but first let us explain the relation between 
vortices and magnetic monopoles, thereby making precise the definition of the orientation of a vortex. 



On the lattice, monopoles are located by so-called Abelian projection. First, one fixes the links up to a residual 
[/(l)-symmetry, which corresponds to an Abelian gauge theory. For example, we can use the maximal Abelian 
gauge [32 to rotate the color vector of the links as much as possible in, say, the cr 3 -direction by maximizing 

R=^2tT[U fl (x)<T 3 Ul(x)cr s }. (2) 

Afterwards, the link variables are replaced by their diagonal part 

Un(x) = anl + id ■ a —¥ u u (x) = — \a a l + ia-icr-i] = I ^ ^ '' r .„ , >-, I . (3) 
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Consider a 3-dimensional cube on the lattice. Normally, the total magnetic flux out of the cube vanishes due to 
div B = 0. If the cube contains a monopole, there would be a nonzero net flux of the (unphysical) Dirac string, which 
compensates the physical monopole flux. The magnetic charge inside the cube is defined by |33j 
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9 p is the sum of the angles 0^ around one plaquette p; ^2 p is always zero as a consequence of the Bianchi identity. 
p is 9p + 2irk with k such that P falls into the range [— 7r, n]. If the absolute value of all plaquette angles is smaller 
than 7r, p = Op and m — 0, as usual. Plaquettes which are greater than 7r are pierced by the Dirac string, at the end 
of which a monopole is located. P discards the flux caused by the Dirac string and a value m ^ results. 

Now let us return to center vortices. Traversing a thick vortex sheet, the link variables change gradually, building 
up to a center element. Pictured in group space of SU(2) ~ S 3 , we travel along a path from unity to the antipode 
— 1. After Abelian projection, only one possible direction for the path remains and this path will go either in +03 
or in —0-3 direction. We can use this sign to allocate an orientation to every patch of the vortex surface. The 
sign also corresponds to the direction of the center flux bundled within the vortex, which is quantized in units ±77. 
Where regions of opposite orientation touch, the flux jumps by 27r, indicating the presence of a sink or source, i.e., a 
magnetic monopole carrying a quantized magnetic charge. To summarize, a center vortex can be imagined as a chain 
of magnetic monopoles, whose flux is bundled within the vortex surface. The monopole worldlines divide the vortex 
sheet into patches of different orientations. 

A. Plane vortices 

In SU{2) lattice gauge theory, plane vortices parallel to two of the coordinate axes are defined by links varying in a 
U(l) subgroup of SU{2). This U(l) subgroup is generated by one of the Pauli matrices o*, i.e., U^ — exp(i(j)ai). The 
direction of the flux and the orientation of the vortices are determined by the gradient of the angle <j), which we choose 
as a piecewise linear function of the coordinate perpendicular to the vortex. The explicit functions for </> are given in 
Eq. (2.1) of |Mj. Upon traversing a vortex sheet, the angle </> increases or decreases by 7r within a finite thickness of 
the vortex. Center projection leads to a (thin) P-vortex at half the thickness [18]. If we consider these thick, planar 
vortices intersecting orthogonally, see Fig. [Ik, each intersection carries a topological charge with modulus \Q\ = 1/2, 
whose sign depends on the relative orientation of the vortex fluxes [5S] . 

B. The Colorful Spherical Vortex 

The nonorientable spherical vortex of radius R and thickness A was introduced in 26 and analyzed in more detail 
in [27] and [29]. It is constructed with the following SU(2) links: 

u{xn= Up(ia(\f-f \)fi-a) « = *,/* = 4 w . th a{f)ss f-K (5) 
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where r is the spatial part of x v and the profile function a is either one of a+ , a_ , which are defined by 
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FIG. 1. a) Two plane vortex pairs intersecting in four space-time points, b) The color structure of the vortex surface from the 
hedgehog configuration leads to c) monopole lines after Abelian projection in the maximal Abelian gauge. 

This means that all links are equal to 1 except for the t-links in a single time slice at fixed t — i^ . The phase changes 
from to Ti from inside to outside for a+(r) (or vice versa for a_(r)). The graph of a~(r) is plotted in Fig. 2 in 26J, 
giving a hedgehog-like configuration, since the color vector n points in the "radial" direction r/r at the vortex radius 
R, see Fig. [T]d. The check that this configuration is a vortex is done with maximal center gauge fixing and center 
projection and results in a P-vortex forming a lattice representation of a 3-sphere of radius R at time slice ti. The 
color structure of the thick vortex leads to a monopole loop on a great circle of the P-vortex after fixing to maximal 
Abelian gauge and Abelian projection. The direction of the loop depends on the U{1) subgroup chosen as Abelian 
degrees of freedom. For the subgroup defined by the Pauli matrices <7i,CT2 or <7 3 the monopole loops are in the yz-, 
zx- and xy-plane, respectively. This is indicated schematically in Fig. Hk with three colors. 

The hedgehog-like structure is crucial for our analysis. The t-links of the spherical vortices fix the holonomy of the 
time-like loops, defining a map Ut(x,t — ti) from the jcyz-hyperplane at t = ti to SU(2). Because of the periodic 
boundary conditions, the time slice has the topology of a 3-torus. But, actually, we can identify all points in the 
exterior of the 3 dimensional sphere since the links there are trivial. Thus the topology of the time slice is R 3 U {oo} 
which is homeomorphic to S 3 . A map S 3 — > SU(2) is characterized by a winding number 
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resulting in N = — 1 for positive and N = +1 for negative spherical vortices. Obviously such windings, given by 
the holonomy of the time- like loops of the spherical vortex, influence the Atiyah-Singer index theorem J5l-f8j. 136J. I37J . 
giving a topological charge Q = — 1 for positive and Q = +1 for negative spherical vortices (antivortices) . Hence, 
spherical vortices attract Dirac zero modes similar to instantons. In |29j we showed that the spherical vortex is in 
fact a vacuum-to-vacuum transition in the time direction which can even be regularized (smoothed out in the time 
direction) to give the correct topological charge also from gluonic definitions (see also [38] for more details). 

III. DIRAC EIGENMODES 

According to the Banks-Casher analysis [39] , chiral symmetry breaking (xSB) is necessarily associated with a finite 
density of near-zero eigenmodes of the chiral-invariant Dirac operator, resulting in a finite chiral condensate, the 
order parameter of chiral symmetry breaking. We compute the lowest-lying chiral eigenvectors xr,l and eigenvalues 
|A| <G [0, 1] of the overlap Dirac operator 
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with the kernel Wilson Dirac operator Dw 40J . Henceforth we will simply write A instead of |A| and assume it 
to be the absolute value of the two complex conjugate eigenvalues of D ov if A ^ 0, 1. We will however distinguish 
between right- and left-handed zero modes of D ov . When speaking of an eigenmode, we always mean both of the 
eigenvectors ip± belonging to one value of A ^ 0, which have the same scalar and chiral densities. For convenience, 



we will number the eigenf unctions in ascending order of the eigenvalues. #0+ denotes the right-handed zero mode, 
#0- the left-handed zero mode. #1 labels the lowest nonzero mode, #2 the second lowest etc. Their eigenvalues are 
referred to as A#0+, A#0- A#l, etc., and their densities as p#0+, etc. Finally, we will name near-zero modes which 
emerge from would-be zero modes also by #0 for good reasons, which will be discussed later. 
The corresponding eigenvectors are given by 

ip± = -7s(XR±iXL)- (8) 
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They have scalar and chiral densities 
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The chiral density p$ is important to assess the local chirality properties, in particular to test the notion that the 
near-zero modes arise from the splitting of exact zero modes localized at lumps of topological charge. 

In a gauge field with topological charge Q 7^ 0, D ov has \Q\ exact zero modes with chirality — signQ, and an equal 
number of eigenvectors of opposite chirality and eigenvalue 1 (doubler modes). This is required in order that Ti'75 = 0. 
The topological charge is 



Q = Tr(7 5 D ov ) = n. 
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(13) 



as for any Ginsparg- Wilson operator. 

The plot titles of the density plots (see e.g. Fig. [2]) give the x- and y-coordinates of the shown zi-slice, the chirality 
(i.e., "chi=0" means we plot ps, "chi=l" would be p+ and "chi=-l" p_), the numbers of plotted modes ("n=l-l" 
means we plot p#l, "n=l-2" would be p#l + p#2) and the maximal density in the plotted area ("max=..." - some 
maxima are cut off in order to resolve other substructures) . For instanton and spherical vortex configurations we use 



12 J x 24-lattices and always plot xi-slices at y 



6 since they are symmetric in spatial directions around their 



centers at x — y 



6.5. For plane vortex configurations we use 16 -lattices; the plotted slices vary. 



A. Eigenvalue Spectrum of the free Dirac operator 



In this section we will analytically calculate the eigenvalues of the massless free Dirac operator. These eigenvalues 
and their multiplicity will come handy later when we compare the eigenvalues for different gauge field configurations 
to the eigenvalues for the free case. The free continuum Dirac operator for massless fermions is given by D = 7^9^, 
where 7^ are the Euclidean Dirac matrices. We want to solve the Dirac equation D a p ipp{x) — A ip a (x). Using the 
well know ansatz of plane wave functions ip a (x) = u a exp(ip M a; /i ), one gets the eigenvalues of the free continuum 
operator A = ±i^/pjj)^. Identifying the eigenvalue A with the fermion mass M makes clear, that this is simply the 
relativistic energy momentum relation in Euclidean space given by M 2 = E 2 + p 2 , where we have denoted p 2 by E 2 
and piPi by p 2 . Note that both the eigenvalue for the plus and the minus sign have a degeneracy of two, which gives 
in total four eigenvalues for the 4x4 Dirac matrix (4x4 Dirac indices) . Additional degeneracy comes from the color 
indices. The free Dirac operator is color blind. Therefore, the degeneracy of the different eigenvalues is multiplied by 
n co i, where n co [ is the range of the color indices. For the overlap Dirac operator the eigenvalues can also be calculated 
analytically for the free case using the lattice ansatz if) a (ri) — u a exp(z a p^n^), with the lattice sites n^ and yield 
the continuum solutions with corrections of orders of a 2 : 



n ( n \- 1 1 , »7m sin(p M ) - MjpJ 
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with M(p^) = M — J2n (1 — C0S (P/j))- Here a — 1 and M = 1 in the free case. This gives {H 
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with expected multiplicities (H 2 v (p^) is a diagonal 4x4 matrix). For p p << 1 we get 

A 2 = Jgf(l + 0( W )). (16) 

With the normalization of Eq. ufj, which leads to A 2 £ [0, 1], a wave function renormalization is needed to convert 
to the usual continuum normalization. This compensates the factor j^p, see e.g. Eq. (6) in (41) . giving the usual 
free Dirac eigenvalues in the continuum limit. Note that the plane wave ansatz is periodic, not only in x^ or lattice 

indices n^, but also in p^. This means, that p^ and p^ H - (with 2„ £ Z) correspond to the same eigenfunction. 

Therefore, in order to get the correct multiplicities for the eigenvalues, we have to restrict the range of p^ to 

- - < P„ < ~ ■ (17) 

a a 

For the usual periodic boundary conditions in spatial directions and periodic or antiperiodic boundary conditions in 
temporal direction, the allowed values for p^ are 
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where N sp is the spatial and N t the temporal extent of the lattice. The total multiplicity of an eigenvalue A is given 
by 

n mu it{y> = 2 n col rip(X), 

where n p (X) is the number of different p M corresponding to a particular A. The factor two comes from the Dirac indices 
as discussed above. Let us now have a quick look at the multiplicity of the eigenvalues for one particular lattice size, 
i.e., N sp — 12 and Nt = 24, which we will use later for our lattice configurations. We assume antiperiodic boundary 
conditions in temporal direction, a — 1 and a SU(2) gauge field theory (i.e., n co i = 2). The momentum vectors 
corresponding to the lowest eigenvalues are then given by pi = and p 4 = ±7r/24. Therefore, we have n p (A#l) = 2 
for the lowest eigenvalues A#l. That means, we get a multiplicity n mM / t (A#l) = 2 n co i n p (A#l) = 8. For the second 
lowest eigenvalues A#2 we have pi = and p^ = ±37r/24 = ±7r/8 and therefore another degeneracy of eight. Then we 
get pi — ±7r/6 (with i £ {1,2,3}) and p^ — ±7r/24, i.e., n p (A#3) = 12. This gives a multiplicity n„ m ; t (A#3) = 48. 

Fig. [2] shows the chiral density of free overlap eigenmodes obtained numerically using the MILC code. The modes 
are found with the Ritz functional algorithm [42, 43] with random start and for degenerate eigenvalues the eigenmodes 
span a randomly oriented basis in the degenerate subspace. Therefore the numerical modes presented in Fig. [2] are 
linear combinations of plane waves with ip^ and show plane wave oscillations of 2p^ in the chiral density. The 
first eight degenerate modes consist of plane waves with p 4 = ±7r/24, hence there is one sine (cosine) oscillation in 
time direction, the next eight have p^ = ±37r/24, i.e., three oscillations in the time direction. The oscillations of 
Xr and Xl are separated by half an oscillation length, i.e., the maxima of p + correspond to minima of p- and vice 
versa. Accordingly, the scalar density p(x^) = \{Xr{ x h)xr{ x n) + Xl^/Jxl^^)) = 1/iVy is constant (Ny ■ ■ ■ lattice 
volume) . 

B. Zero Modes and Instantons 

An instanton field gives rise to an exact zero mode (the precise expression is given in e.g. [44] ) . Its probability 
density 

2 R 2 

p{x) = ^( X - Xo r + Ri (18) 

is localized at the instanton core Xo, with a half-radius of the instanton parameter R. Since the zero mode is exactly 
chiral, its chiral density 'i/^75'0 equals ±p. 

In the instanton liquid model the near-zero modes originate from the overlap of the would-be zero modes carried 
by individual instantons and anti-instantons. If the overlap is not too large, one expects that the resulting near-zero 
modes still exhibit definite chirality locally. This picture predicts characteristic properties of the low-lying modes [45] : 

• Their probability density should be clearly peaked, indicating the location of instantons. 



y=6, z=6, chi=0, n=l-l, max=0.0000 149759 



y=6, z=6, chi=0, n=7-7, max=0.0000206666 



y=6, z=6, chi=-l, n=9-9, max=0.0000174434 




FIG. 2. Chiral density of the low-lying eigenmodes of the free overlap Dirac operator: ps#l (left), ps#7 (center) ps#9 (right). 
The modes clearly show the plane wave behavior with oscillations of 2p M (see text). 

• The local chirality at the peaks should match the sign of the topological charge and the size of the chiral lump 
should be correlated to the extension of the topological structure. 

• As an instanton and an anti-instanton approach each other, the eigenvalues should be shifted further away from 
zero and the localization and local chirality properties should fade. 

Numerical evidence supporting these assumptions about the local chirality structure of the low-lying Dirac modes 
are presented in e.g. [3B] and [37]. Here we want to analyze these issues for center vortices. 



Zero Modes and Center Vortices 



Reinhardt et al. [35] analytically calculate the exact zero modes of the Dirac operator in the background of plane 
vortices, both nonintersecting and perpendicularly intersecting ones, for an Abelian £7(1) gauge field on R 2 ,R 4 ,T 2 
and T 4 . 

The first vortex field of interest consists of two (anti)parallel fluxes on a two-torus. Because of a two-dimensional 
translation symmetry we can identify this with the four-dimensional configuration of a single vortex pair presented in 
Sec. |II A| The second example contains four flat vortices on a four-torus, which intersect orthogonally in four points. 
This corresponds to the configuration shown in Fig. [Tk. 

We want to discuss the boundary conditions used in [48j . The four-dimensional problem can be reduced to two 
two-dimensional ones because of the translation invariance parallel to the vortex sheets. It therefore suffices to deal 
with the T 2 case. The boundary conditions for a gauge field on T 2 are specified by the two transition functions 
£l x ,tt y . In a suitable gauge, it is possible to set ft x = 1. The co-cycle condition then turns into a periodicity condition 
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(19) 



for £l y (x), which lives only at the boundary of the rectangle which represents the torus. Therefore Q y defines a 
function S 1 — > S 1 belonging to a class of the homotopy group 7r 1 (S' 1 ) ~ 7L. Its winding number n determines the 
magnetic flux through the xy-plane by 



$:=</> A^dx^ — 2-kti . 



(20) 



To prove this, we simply perform the line integral over the boundary of the rectangle resulting from cutting up the 
torus. We set £l y (x) = exp{i\(x)}. Since Q y is periodic, x( x + L x ) = x( x ) + 27m. Then 

f 1 ! rL y rL x rL y 

$= / A x {x,0)dx+ / A y (L x ,y)dy- / A x (x, L y )dx - / A y (0,y)dy 
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(21) 
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FIG. 3. a) The lowest eigenvalues for the spherical vortex with R = d = 3a compared to the eigenvalues of the free Dirac 
operator (red crosses). Eigenvalues are calculated on a 12 3 x 6 lattice with relative lattice constant a — 2 (magenta boxes), on 
a 16 3 x 8 lattice with a = § (blue stars) and on a 24 3 x 12 lattice with a = 1 (green crosses). The multiplicity of A#l is eight, 
of A#2 48. b) Spatial distribution of the zero mode, mainly located at the vortex core. 

Consequently, the flux on T 2 is quantized. Note that this would not hold for SU(2) since iri(S 3 ) — 1. To create two 
vortices which carry a flux of +ir each, Reinhardt et al. use Q y — exp{27rix}. Returning finally to the fermion field, 
it obeys boundary conditions consistent with gauge invariance, 



ip(x, y + L y ) = n y (x)tp(x, y) . 



(22) 



In [34] we analyzed the zero modes for plane vortex configurations and found good agreement with the results 
obtained by Reinhardt et al. |48) . The discrepancies discussed in our work originate in the finite thickness of our 
vortex configurations because of finite lattice sizes. 

We conclude this section with a quick discussion of the eigenvalues for the spherical vortex. As already mentioned, 
there always occurs exactly one zero mode, a positive chirality one for the spherical vortex and a negative chirality 
zero mode for the antivortex. The lowest nonzero modes can be seen as some sort of modified eigenmodes of the free 
Dirac operator. This point of view is motivated by the results presented in Fig. [3k. In this diagram, one can see the 
eigenvalues for configurations with vortices of the same size in lattice units but different spatial N sp and temporal 
N t lattice extents. In other words, the vortex gets smaller while the lattice gets finer. One sees from the figures, 
that the nonzero eigenvalues converge to eigenvalues of the free Dirac operator as the vortex gets smaller and smaller. 
For the investigated eigenvalues, there seems to be a one to one correspondence between the eigenvalues of the free 
Dirac operator and the nonzero eigenvalues of the Dirac operator for the vortex configuration. The zero mode is 
not a lowered nonzero mode, it occurs in addition to the low lying modes. Because the total number of eigenmodes 
only depends on the lattice size, the number of complex eigenvalues is lowered by two (because of the zero and the 
doubler mode) for the vortex configuration in comparison to the free case. Note that we still can have a one to one 
correspondence between all the complex eigenvalues for the vortex and the free case. One can see this by remembering 
that the correspondence is established in the limit of an infinitely large lattice. 



IV. INTERACTIONS BETWEEN TOPOLOGICAL OBJECTS 



We want to discuss the Dirac equation for a gauge field An = A\^ + A2 ii consisting of two fields A\^ L and y^ 
which are separated in Euclidean space and have nonvanishing topological charge Qi. Therefore, the Dirac operators 
D\ for Aifj, alone and D 2 for A 2l i alone would have at least one zero mode. In the following, the zero modes of D\ 
and D 2 will be called would-be zero modes. Lets discuss the case in which Ai^ has Q\ = 1 and A 2il has Q 2 = — 1. 
For simplicity we assume that Di has only one left handed zero mode | tpi) and D 2 only one right handed zero mode 
| ip 2 ). Clearly, this two would-be zero modes are orthogonal, they can be part of an orthogonal basis. Let us now have 
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a look at the Dirac operator in this orthogonal basis. In particular we are interested in the upper left 2x2 block 

(23) 
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of the Dirac matrix. The continuum Dirac operator D is given by 

D = 7 M [dp + i^( x )) = 7m (#m + ^imW + ^^(aO) 
= Di+ jf 1 iA 2 ^(x) = D 2 + 7m i-^iM^) ■ 



(24) 



The first element of ( 23 1 therefore evaluates to 

(^ | D | 0i) = (0i I A I 0i) + (01 I 7m WjmW I 00 = • 
One can see that (^i | 7 M iA2n(x) \ ipi) vanishes from 

(tpi | 7^ L4 2M (a;) I 01 ) = ("01 I 7m iA^O) 75 2 I 0i) 

= - ("01 I 75 7m i-^m^) 75 | 0i) = - (01 I 7m [ -A 2 n(x) | 0i) , . 
In the same way one can prove 

(V* I -D | V»> = . 
Let us now calculate the off-diagonal terms. The first off-diagonal term evaluates to 

(fh\D\ V> 2 ) = (0i | D 2 | V 2 ) + (0i | 7m Ui/i(*) I 02> = + c = c. (28) 

Here c stands for the overlap integral (tpi | 7^ L4 lAt (x) \ 1P2). In general, c will be large for eigenmodes that overlap 
a lot, and small for eigenmodes that overlap only a little. Note that it is crucial that | ipi) and | 1P2) have different 



(25) 



(26) 



(27) 



chirality. Otherwise, by the same argument as used in (26), the overlap integral would have to vanish. The second 
off-diagonal term of the upper-left 2x2 block is 



(02 I D I 0i) = ((0i I £>* I 2 })* = (-{ih\D\ 2 ))* = -c* . 
Here it was used that the continuum Dirac Operator D is antihermitian, i.e., D^ = —D. 



(29) 



Combining (25), (27), (28) and (29), the upper left 2x2 block of the Dirac matrix reads 



c 
-c* 



(30) 



This 2x2 block can easily be diagonalized. The eigenvalues Ai,2 and the normalized eigenvectors tl>[ 2 of (30) are 



Ai, 2 = ±h 



= ±i|c| 



01,: 



1 

7^ 



±i 



1 



(31) 



This means that the interaction transforms the two (would-be) zero modes into two near-zero modes. Those near-zero 
modes also occur additionally to the free Dirac eigenmodes and therefore we still enumerate them like zero modes by 
#0. The strength of the interaction quantified by the overlap integral c determines the size of the near-zero eigenvalue. 
Note that the new near-zero modes consist in equal parts of the would-be zero modes | ipi) and | 1^2) ■ Therefore the 
scalar and chiral densities of the near-zero modes are simply an average of the densities of the would-be zero modes. 
So far we have ignored everything except the upper left 2x2 block of the Dirac matrix. To get exact eigenmodes 
we clearly have to diagonalize the whole Dirac matrix and not only this 2x2 block. However, (31) represents a 
legitimate approximation to the exact eigenvalues and eigenmodes if the elements of the form 

(0j I D I ip!.?) = (tpj | 7m iAimW I 0i,2) with j > 2 

are a lot smaller than the overlap integral c = (V>i | 7 M iA\^(x) \ ^2)- Usually this will be the case, because | ip\) 
is localized at A\^, and | ij)j) with j > 2 is not. Let us also have a quick look at what happens when we have two 
would-be zero modes with the same chirality. In this case also the off-diagonal terms of the upper left 2x2 matrix 
vanish and the would-be zero modes are actual zero modes. 

Note that the mechanism discussed in this section is the basis for the instanton liquid model of spontaneous 
chiral symmetry breaking (see [H] for a detailed review). In the instanton liquid model the QCD- vacuum consist 
of an ensemble of instantons and anti-instantons whose would-be zero modes split into near-zero modes because of 
interactions. Therefore, one gets a nonvanishing eigenmode density around zero, which gives via the Banks-Casher 
relation a finite chiral condensate and broken chiral symmetry. Clearly, one can construct such a model also with 
other topological objects, as will be shown in the next section for center vortices. 
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FIG. 4. The lowest overlap eigenvalues for instanton and spherical vortex configurations compared to the eigenvalues of the 
free (overlap) Dirac operator. 



V. VORTICES AND X SB 



In Sec. |TT] we discussed the various possibilities of vortices to create topological charge, i.e., via writhing and 
intersection points as well as through their color structure. In previous works we presented results on the attraction 
of zero modes for flat [31] and spherical [2l)lI29| vortex configurations. Here we present some results on how vortices 
form near-zero modes from would-be zero modes through interactions. 



A. Spherical Vortices and Instantons 



We start with spherical vortices and show that their effects on fermions is pretty much the same as those of 
instantons. We have shown in Fig. [3] that they attract a zero mode and its scalar density peaks at the vortex surface. 
We interpreted the nonzero modes as eigenmodes of the free Dirac operator, which are shifted slightly because of their 
interaction with the nontrivial gauge field content. In Fig. [4^, we see that a single instanton has nearly exactly the 
same eigenvalues as a single spherical vortex. In Fig.[5]we show that even the chiral densities of the lowest eigenmodes 
distribute similarly, except for the fact that the response of the fermions to the spherical vortex is squeezed in the 
time direction, since the vortex is localized in a single time slice (t — 5) . Another interesting issue is that the nonzero 
eigenmodes show nice plane wave oscillations, like the free eigenmodes in Fig. |2j mode #8 however shows some 
similarity to the zero mode, its eigenvalue is also clearly enhanced compared to the free spectrum. This is only a side 
remark however, as we are not sure how to interpret this and it does not seem to be important for the creation of 
near-zero modes since we observe the same effect for instantons. 

We further plot the spectra of instanton-anti-instanton, spherical vortex-antivortex and instanton-antivortex pairs 
in Fig. Hk. We again see nearly exactly the same eigenvalues for instanton or spherical vortex pairs, but now we 
get instead of two would-be zero modes a near-zero mode for each pair. The chiral density plots in Fig. [6] for the 
instanton-anti-instanton pair and Fig. [7] for the spherical vortex-antivortex pair show, besides the similar densities, 
that the near-zero mode is a result of two chiral parts corresponding to the two constituents of the pairs. The nonzero 
modes can again be identified with the free overlap modes with the same side remark for mode #8. 

Finally in Fig. [4b we plot the eigenvalues of two (anti-)instantons and two spherical (anti)vortices giving topological 
charge Q = 2 (Q = —2) and therefore two zero modes, two vortex-antivortex pairs with two near-zero modes and a 
configuration with two vortices and an antivortex (i.e., a single vortex plus one vortex-antivortex pair) giving one zero 
mode (Q = 1) and one near-zero mode. The chiral densities for the last configuration in Fig.^lstiow that the zero mode 
peaks at both spherical vortices, the near-zero mode again consists of two chiral parts from the (second/would-be) 
zero mode of the two vortices and the (would-be) zero mode of the antivortex. Now the modes #7 and #8 clearly 
deviate from the free eigenvalues showing similar densities as the zero and the near-zero mode respectively. The 
results clearly show that we may draw the same conclusions for spherical vortices as for instantons concerning the 
creation of near-zero modes. 
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FIG. 5. Chiral densities of overlap eigenmodes: a) zero mode (left), first (center), ninth (right) and b) eighth (p$ left, p+ center 
and p- right) nonzero modes for an instanton; c) and d) the same as a) and b) but for a spherical vortex. 
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FIG. 6. Chiral densities (ps left, p_ center and p+ right column) of the a) lowest (near-zero), b) second-lowest (nonzero) and c) 
eighth (nonzero) eigenmode of the overlap Dirac operator for an instanton-anti-instanton pair, d) p$ of the sixth (left), seventh 
(center) and ninth (right) eigenmode. 
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FIG. 7. Same as Fig. ml but for a spherical vortex-antivortex pair. Chiral densities (ps left, p+ center and p_ right column) 
of the a) lowest (near-zero), b) second-lowest (nonzero) and c) eighth (nonzero) eigenmode. d) ps of the sixth (left), seventh 
(center) and ninth (right) eigenmode. 
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FIG. 8. Chiral densities of overlap eigenmodes for spherical vortices at t — 4 and t = 20 and an antivortex at t = 12: a) zero 
mode /9_#0— (left), ps#l (center) and ps#5 (right), b) near-zero (second-lowest) mode #0, c) #7 (ss #0—) and d) #8 (« #0) 
with ps (left), p+ (center) and p_ (right). 
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B. Planar vortex configurations 



From Sec. Ill C and Rcinhardt et al. [48] we recall that plane vortices attract zero modes according to their magnetic 



flux. For a single pair of parallel vortex slices with the same flux direction we may get two zero modes, peaking at 
the two vortex slices with the same (local) chirality, see Fig. 10 1. But in four dimensions they can be removed by 
antiperiodic boundary conditions in perpendicular directions of the vortex flux. However, if we intersect two such 
vortex pairs orthogonally, we get two real zero modes, which we analyzed in |34j . Now, for pairs of (anti)parallel 
vortex slices with opposite flux directions the situation is similar, except that we only get near-zero modes, it is 
not clear though how to split them into would-be zero modes. For a single vortex pair we get four near-zero modes 
peaking at the two vortex slices with opposite (local) chirality, see Fig. 10 3. The nonzero modes again show plane 
wave oscillations, see Fig. |10b . If we intersect antiparallel vortex pairs orthogonally we get two intersection points with 
topological charge Q = +1/2 and two intersection points with topological charge Q = —1/2. For such a configuration 
we get four real near-zero modes, with local chirality peaks at the intersection points, according to their topological 
charge contribution, see Fig. [9b. Again, the near-zero modes cannot be split into would-be zero modes which could 
be assigned to single vortex intersections, but the modes rather peak at all intersection points. Hence, the mechanism 
of Sec. |VA"] does not apply in this case, since there are no localized lumps of topological charge Q = ±1. Nevertheless 
the vortices attract chiral near-zero modes via their intersections with topological charge Q — ±1/2. We expect that 
vortex intersections, writhing points and even color structure contributions of vortices to topological charge are able 
to create a finite density of near-zero modes and break chiral symmetry via the Banks-Casher relation. 
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FIG. 9. a) The lowest overlap eigenvalues for plane vortex configurations compared to the eigenvalues of the free (overlap) 
Dirac operator (red crosses), b) Chiral density ps#0 in the intersection plane of four near-zero modes of crossing flat vortex 
pairs with opposite flux direction (Q = 0). 
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c) 

FIG. 10. Chiral densities of the low-lying eigenmodes of the overlap Dirac operator: a) two zero modes of one plane vortex 
pair, i.e., two vortex slices with the same flux direction, b) four near-zero modes of one plane vortex pair, i.e., two vortex slices 
with opposite flux direction and c) first nonzero mode of the crossing flat vortex pairs in Fig. [9b). 
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C. Asqtad Staggered Modes 



For completeness we shortly discuss the asqtad staggered eigenmodes for the presented configurations. In principle 
the same conclusions as for overlap modes apply if we consider the double degeneracy of asqtad staggered modes 
due to charge conjugation. Hence we have two (would-be) zero modes for Q = ±1 and there are four times as many 
near-zero modes since two would-be zero modes do result in two near-zero modes instead of one for overlap modes. We 
further plot the chirality (i^j^t/j) of the asqtad staggered eigenmodes, where 75 corresponds to a displacement along 
the diagonal of a hypercube. To ensure gauge invariance the product includes gauge field multiplications along all 
shortest paths connecting opposite corners of the hypercube. Staggered fermions do not have exact zero modes, but a 
separation between would-be zero modes and nonchiral modes is observed for improved staggered quark actions 49J. 
The plots show that the would-be and near-zero modes have enhanced chirality compared to nonzero modes and we 
even observe the local chiral density properties for the near-zero modes which we discussed for the overlap modes. 
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FIG. 11. The lowest asqtad staggered eigenvalues for instanton and spherical vortex configurations compared to the eigenvalues 
of the free Dirac operator, b) Chirality of the corresponding eigenmodes. 
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FIG. 12. a) The lowest asqtad staggered eigenvalues for plane vortex configurations compared to the eigenvalues of the free 
Dirac operator, b) Chirality of the corresponding eigenmodes. 
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FIG. 13. Chiral densities of the low-lying eigenmodes of the asqtad staggered Dirac operator: a) four near-zero and b) non-zero 
modes of a spherical vortex-antivortex pair, and c) 16 near-zero modes of the crossing flat vortex pairs in Fig. Md). 



VI. CONCLUSIONS 



The instanton liquid model provides a physical picture of chiral symmetry breaking via the idea of quarks hopping 
between random instantons and anti-instantons, changing their helicity each time. This process can be described by 
quarks propagating between quark-instanton vertices. As fermions do not seem to make much of a difference between 
instantons and spherical vortices this picture can be extended to colorful spherical center vortices. But in the vortex 
picture the model of chiral symmetry breaking can be formulated even more generally, as we have shown that various 
shapes of vortices attract (would-be) zero modes which contribute via interactions to a finite density of near-zero 
modes with local chiral properties, i.e., local chirality peaks at corresponding topological charge contributions. The 
simple picture of localized would-be zero modes from the instanton liquid model or spherical vortex configurations as 
discussed in Sec. IV and V A does not apply directly to general vortex structures as there are not only topological 
charge contributions of Q = ±1. In Monte Carlo configurations we do not, of course, find perfectly flat or spherical 
vortices, as one does not find perfect instantons. The general picture of topological charge from vortex intersections, 
writhing points and even color structure contributions or instantons can provide a general picture of x^B: any source 
of topological charge can attract (would-be) zero modes and produce a finite density of near-zero modes leading to 
chiral symmetry breaking via the Banks-Casher relation. Here one also has to ask what could be the dynamical 
explanation of x^B. We can try the conjecture that only a combination of color electric and magnetic fields leads 
to xSB, electric fields accelerating color charges and magnetic fields trying permanently to reverse the momentum 
directions on spiral shaped paths. Such reversals of momentum keeping the spin of the particles should especially 
happen for very slowly moving color charges. Alternatively we could argue that magnetic color charges are able to 
flip the spin of slow quarks, i.e. when they interact long enough with the vortex structures. 

Finally, it seems that vortices not only confine quarks into bound states but also change their helicity in analogy 
to the instanton liquid model. In accordance with Casher's argument, a force strong enough to confine quarks is 
also generally expected to break chiral symmetry [50] , we therefore conclude that the center vortex model of quark 
confinement may indeed be capable of describing chiral symmetry breaking. We cannot give a conclusive answer 
to the question of a dynamical explanation for the mechanism of x^B and only speculate on the importance of our 
results for Monte Carlo configurations, since there vortices are neither perfectly flat nor spherical, as there are no 
perfect instantons either. But the importance of the long-range nature of low-dimensional topological structures for 
the understanding of the mechanism of xSB in QCD was underlined by various results of different groups [5TM57J . 
and agrees well with a vortex picture of %SB. 
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